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Abstract 

Distribution of water concentration on the Moon is currently an area of active research. Recent 
studies suggest the presence of ice particles, and perhaps even ice blocks and ice-cemented regolith on the 
Moon. Thermal extraction of the in-situ water is an attractive means of satisfying water requirements for a 
lunar mission. In this paper, a model is presented to analyze the processes occurring during the heat-up of 
icy regolith and extraction of the evolved water vapor. The wet regolith is assumed to be present in an 
initially evacuated and sealed cell which is subsequently heated. The first step of the analysis involves 
calculating the gradual increase of vapor pressure in the closed cell as the temperature is raised. Then, in 
the second step, the cell is evacuated to low pressure (e.g., vacuum), allowing the water vapor to leave the 
cell and be captured. The parameters affecting water vapor pressure build-up and evacuation for the 
purpose of extracting water from lunar regolith are discussed in the paper. Some comparisons with 
available experimental measurements are also made. 

Nomenclature 

A number of molecules per unit mass of water (kg' 1 ) 

AH enthalpy of vaporization (Jkg' 1 ) 

J s water vapor flux (kgm'V 1 ) 

M molecular species in bed with attached water of hydration 

P pressure (Pa) 

P e vapor pressure in reactor from evaporation due to heating (Pa) 

P r regulated reactor pressure (Pa) 

P s vapor pressure at the condensed-phase surface (Pa) 

P 0 vapor pressure in reactor prior to heating (Pa) 

R g gas constant for water (Jkg' 1 K' 1 ) 

S surface area for vaporization (m 2 ) 

T temperature (K) 

T r reactor temperature (K) 

T s vapor temperature at the condensed-phase surface (K) 

V free volume in reactor (m 3 ) 

f( V)dv Maxwellian velocity distribution function (m" 3 ) 

k mass transfer coefficient (ms' 1 ) 

m mass of condensed-phase water in bed (kg) 

m 0 mass of condensed-phase water in bed prior to heating (kg) 

772, mass of water vapor vented during the initial evacuation period (kg) 

n s number density of water vapor molecules near the condensed-phase surface (m' 3 ) 

t time (s) 
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t c time constant for evaporation (s) 

v velocity (ms' 1 ) 

k e evaporation coefficient 

p 0 density of water vapor in the reactor prior to heating (kgm' 3 ) 

p, density of water vapor at the condensed phase surface (kgm' 3 ) 

p,. density of water vapor in the reactor (kgm' 3 ) 

Introduction 

A major objective of lunar in-situ resource utilization (1SRU) (Refs. 1 to 3) is production of water and 
oxygen from lunar regolith. To this end, the distribution of water concentration on the Moon continues to 
be an area of active research. Hydrogen concentrations, at least some of which were presumed to indicate 
the presence of water, were observed near the lunar poles using remote neutron scattering (Ref. 4). More 
recently, the Lunar Crater Observation and Sensing Satellite (LCROSS) mission (Ref. 5) indicated water 
concentrations in the range of 5.6 ± 2.9 percent by mass at the impact site in the permanently shadowed 
region of the crater Cabeus at the lunar south pole. It is believed that the ice may be present in the form of 
ice particles, ice-cemented regolith and even ice blocks. Away from the poles and in sunlit areas of the 
Moon, spectroscopic measurements have shown hydroxyl and water absoiption at many locations. Since 
any surface water in regions exposed to the sun is unstable and would have sublimated, the observed 
signals are interpreted as arising from water adsorbed at the surface or present within the interstices of the 
lunar regolith. The absorption signals indicate a trapped water content of 1 0 to 1 000 ppm. 

A thermal process where the icy or water-containing regolith is heated to remove the water in the for 
of vapor may be utilized for water extraction. There are differences in the physical processes that occur 
depending upon whether the water is in the form of ice (either as icy particles or regolith-ice blends) or 
bound within the regolith particles. In the case of regolith-ice blends, heat addition primarily creates a 
transformation of the ice to the vapor phase which may be removed. The physical transformation in the 
general case is 


H 2 0(ice) -> H 2 0(g) 

with or without an intermediate liquid phase depending upon the ambient conditions. In the case of bound 
water, the heat addition must first break any bonds between the water molecules and other molecules 
inside the regolith particle (e.g., water may be in the form of water of hydration attached to a different 
molecule), i.e., 


M • (xH 2 0) — » M + xH 2 0 

Subsequently, the water molecules must diffuse out of the particle and then be removed utilizing , for 
example, a carrier gas. The vaporization process for the icy regolith can be carried out at lower 
temperatures (i.e., 100 to 150 °C) than the extraction of bound water which typically requires 
temperatures in excess of 250 °C in order to break the bonds of hydration. 

The NASA project RESOLVE (Regolith & Environment Science and Oxygen & Lunar Volatile 
Extraction) (Ref. 6) is being designed to extract water from both the icy-regolith blends as well as that 
bound within the regolith particles. A batch process is envisioned where a given amount of lunar regolith 
is deposited inside a tubular reactor. The reactor walls are heated to raise the regolith to the desired 
temperature. For the case of the icy regolith, the reactor is vented once there is sufficient vapor pressure 
build-up and the water is captured downstream. For the case of bound water, the regolith is heated as 
before but to a higher temperature as noted above. An inert carrier gas, e.g., argon, is flowed in order to 
remove the water diffusing out of the particles. 

The RESOLVE project will also investigate extraction of water using hydrogen reduction of the 
oxides in the regolith particles. The hydrogen reduction method targets mainly the iron oxide in the 
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regolith. Iron oxide occurs in lunar regolith as ilmenite, particularly in titanium-rich basalts, or may be 
found dispersed in the glassy phase of lunar regolith particles. The overall gas-solid reaction between FeO 
and hydrogen is given by: 


FeO + H 2 <=>Fe + H 2 0 

The formed water may be electrolyzed in an electrolyzer subsystem to produce oxygen and recycle the 
hydrogen. 

While thermal and hydrogen reduction systems are well developed for Earth applications, in the lunar 
context it is important to establish a favorable balance between efficiency, robustness, and equivalent 
system mass in order to conserve resources. The reactor must interface with the other subsystem 
processes such as upstream regolith extraction and downstream electrolysis and phase separation. A need 
exists, therefore, to establish an analytical framework to couple these subsystems together. The building 
blocks of this framework are the individual subsystem models. 

The development of a generalized analytical model of the hydrogen reduction reactor subsystem and 
its validation against experimental data from reduction of JSC-1A, a lunar regolith simulant, has been 
discussed in earlier papers (Ref. 7). The overall objective of this paper is to model water extraction from 
icy regolith in a tubular reactor geometry similar to that envisioned for the RESOLVE project. The 
extraction process is formulated in terms of the following two steps: (i) water vapor evolution and 
pressure build-up in the closed reactor as it is heated, and (ii) venting of the vapor to a secondary reservoir 
maintained at a suitable (lower) pressure. The key parameter appearing in the model is a characteristic 
time (i.e., time constant) for evaporation which depends upon the evaporating surface area, the resistance 
of the regolith bed to vapor transport, and the free volume available to the vapor. 

In the following, the model development is discussed first. Then, calculated results from the model 
are presented. These include the effects of different values of the time constant, reactor temperature and 
vent pressure on the vapor generation and ice depletion rate. Qualitative comparison with initial 
experimental results on the vapor pressure build-up are also provided. 


Model Description 

The water extraction process begins by filling the empty reactor with the desired amount of icy 
regolith. The reactor is then sealed. The walls of the reactor are heated externally. During this heating 
period, the temperature of the regolith increases. Water vapor is given off and accumulates in the free 
volume of the reactor. When the water vapor pressure reaches a pre-determined level or the system 
reaches equilibrium, the vapor is vented into a reservoir where it may be condensed to liquid. The model 
addresses both of these aspects. 


Pressure Rise Due to Vapor Build-Up 

At the surface of the condensed phase (i.e., ice or liquid water) it may be assumed that there is 
thermodynamic equilibrium which yields a unique relation for the vapor pressure in terms of the 
temperature. This is given by the Clausius-Clapeyron relation (Ref. 8) 


dP_ 

P 


A H 

W 1 


dT 


( 1 ) 


The vapor pressure difference at the surface of the condensate and in the free volume of the reactor is the 
driving force for continued evaporation of the condensed phase. The vapor mass transfer from the regolith 
bed to the free volume of the reactor is hindered by the depth of the bed which resists the transport. 
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Thus, apart from the amount of ice present initially in the dirt sample, parameters important for the 
analysis are: 

(a) bed/reactor temperature: increasing the temperature will increase the vapor pressure at the particle 
surface, thereby increasing the driving force for evaporation, 

(b) free volume: increasing the free volume will lead to reduced partial pressure of water vapor in the 
free volume, thereby increasing the driving force for evaporation, and allowing for an increased amount 
of water vapor collection, 

(c) mass transfer coefficient: decreasing the mass transfer coefficient (thereby increasing the 
resistance of the dirt bed to water vapor transport) will decrease the vaporization rate, and 

(d) surface area of the condensed phase: increasing the surface area of the condensed phase will 
increase the area over which vaporization takes place. The overall amount of ice/dirt also contributes to 
the surface area. 

The model considers the regolith bed as a lumped unit which is suitable for well-mixed beds during 
operation. If, for example, a sufficiently deep dirt bed happens to remain stationary during operation then 
the upper layers of the dirt in contact with the free volume may dry out faster, and the analysis of spatial 
variations of the water vapor profile in the bed may be warranted. 

The model is developed by first considering the case where evaporation takes place to vacuum 
(Ref. 9) and then generalizing to a nonzero vapor pressure in the reactor. The number density of the water 
molecules at the condensed phase surface is related to the corresponding vapor density by Avogadro’s 
Law so that 


n s = Ap s (2) 

The single component Maxwellian distribution function yields the number of molecules with velocity v in 
the velocity space dv 


f(v)dv ■■ 


(2nRJ s y /2 


exp 


V 2R g Ts j 


dv 


(3) 


where v represents the component of the velocity in the direction of the outward normal to the 
evaporating surface. This expression may be used to evaluate the number of molecules evaporating from 
the surface. 

Define an evaporation coefficient A e which represents the fraction of molecules that do not return to 
the surface. Then the flux of water vapor from the evaporating surface is given by the half-range integral 


d s — ^ePs f 


(2nR T s ) 


1/2 


exp 


2 R g T s j 


dv 


which simplifies to 


d s — AT' 


R g R s 

2 n 


(5) 


The rate of change of condensed-water mass evaporating into vacuum may therefore be written as 


dm 

dt 


= ~kS p. 


(6) 
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where the factor k is given by 


R,,T, 

k = K 1 s 


2k 


( 7 ) 


For the case when the vapor pressure in the reactor is nonzero, the equation for m is generalized by 
modifying Equation (6) to include this effect, i.e., 


~T~ ~ ~kS{Ps -Pr) 
at 


( 8 ) 


where p s and p r , are the densities of the surface vapor and the bulk reactor vapor, respectively, and are 
related to their separate vapor pressures as well as the overall vapor temperature (assumed to be the same 
as the reactor temperature) by the equation of state. Since vaporization is not freely into vacuum, the 
factor k is identified as a mass transfer coefficient which accounts for the resistance resulting from the 
finite reactor pressure and bed depth and has units of velocity. In general, k, or the time constant, t c , 
defined later, is expected to be determined experimentally. 

If m 0 is the initial mass of condensed water in the dirt and the initial density of vapor in the free 
volume is p 0 , the quantity p r may be expressed as 


Pr 


- Po + 


(r»o ~ m) 
V 


( 9 ) 


and the equation for the mass of the condensed phase in the reactor, m, becomes 


dm kS 


dt V 


m = -kS\ 



m o 

V 


J) 


( 10 ) 


Once m is determined from the integration of Equation (10), the vapor density in the reactor is obtained 
by Equation (9). 

It is assumed that the mass flux of evaporating water is not large enough relative to the thermal mass 
of the regolith simulant to cool the surface temperature due to the removal of the latent heat of 
vaporization. This assumption eliminates the coupling between heat and mass transport processes. If 
needed, however, an energy equation for the temperature change can be coupled with the vaporization 
rate equation. For the purposes of this paper the temperature-time history of the reactor will be assumed to 
be provided. Note that, in general, the quantities k, S, and V are time-dependent. For example, even for a 
fixed reactor temperature, as the ice coating thins out by evaporation, the resistance to vapor transport 
may change thereby modifying the factor k. At the same time, the surface area for evaporation, S, may 
also change.The free volume may also be altered if the bed void fraction changes due to evaporation of 
the condensed water. Flowever, in order to obtain an estimate of the evaporation rate the factors k, S, and 
V may be assumed constant during the evaporation process. If the reactor temperature is also fixed the 
following solution for m is obtained 


(mp ~ m) 

V 


(Ps — Po ) 


1-exp - 


ks Vi 

— t 

V JJ 


It is clear from the above equation that the parameter 


t 


C 


v_ 

kS 


(11) 


(12) 
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is a characteristic time (i.e. time constant) for the evaporation process. 

Assuming perfect gas behavior, the vapor pressure due to evaporation is given by 


P e ={P,-P*) 


1-exp 


kS A 

1 

V 


yj 


or, equivalently, 


Pe 

(Ps-Po) 


f ( 

1-exp 

v 


v 






J 


(13) 


(14) 


In general, however, the vaporization rate equation must be solved numerically if the temperature changes 
during the evaporation process are to be accounted for. 


Evaporation During Reactor Venting 

The mass of vapor available in the free volume at the end of the heat-up cycle is p x V. This amount is 
generally much less than the water still present in the condensed phase. Thus, vaporization from the 
condensed phase will continue during the vent. It is assumed that the reactor is vented at a constant 
pressure, e.g., using a backpressure regulator. Thus, as vapor leaves the reactor, more condensed water 
evaporates to increase the vapor pressure in the reactor to the regulated value. After the initial vent, this 
may be modeled as a continuous process with the rate of evaporation given by Equation (8). 

It is convenient to cast Equation (8) in terms of the saturation and regulated reactor pressures using 
the equation of state to yield 


dm 

dt 


kS 

R g T r 


( Ps-Pr ) 


(15) 


In terms of the time constant t c and the free volume V the above equation becomes 


dm 

dt 


V 

t c R g T r 


( Ps-Pr ) 


(16) 


It is clear that for vaporization to take place the regulated-reactor pressure must be less than the saturation 
vapor pressure at the reactor temperature. 

Let the mass of water removed when the reactor pressure is first relieved be m,. Then, if the time 
constant and V are assumed to be constant during the vent process, the time for complete evaporation is 
given by 


t 


evap 


(fflp ~ mj) 
dm 


dt 


(17) 


Equations (16) and (17) may be written as 


dm 

dt 


V 


(m 0 - mj)R g T r 


(Ps-Pr) 


where the character A designates a nondimensional quantity with 


(18) 
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- _ 1 
tevap ~ ~dkh 
di 


( 19 ) 


and 


m 

m = 


and 


- t 
t = — 

K 

The quantity (m ir m,)/V is the effective density of the condensed-phase water if it occupied the free 
volume V of the reactor. It is expected to be of the order of 100 kg/m 3 for the conditions to be encountered 
by the RESOLVE experiment. 


Results and Discussion 

Effect of the Time Constant on Evaporation 

Equation (14) indicates that the vapor pressure in the closed reactor cannot exceed the saturation 
vapor pressure at constant temperature. The time constant, t c , is a measure of how fast the vapor pressure 
in the reactor increases with time. A plot of Equation (14) is presented in Figure 1. 

As t c increases, it takes longer for the vapor pressure in the reactor to build up. Equation (12) shows 
that t c is proportional to the free volume V and inversely proportional to the mass transfer coefficient, k, 
and the effective surface area for vaporization, S. Thus, increasing the size of the reactor for a given 
loading of wet regolith will increase the time needed for the vapor pressure to grow to a desired level. The 
mass transfer coefficient, k, is a measure of the ease with which the vapor leaves the regolith bed. As k 
increases the time constant decreases and there is a faster build-up of vapor pressure in the reactor. It may 



Figure 1. — Nondimensional vapor pressure increase in the reactor with 
time at a constant temperature. 
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Figure 2. — Increase in reactor vapor pressure with time during heating 
for two different time constants. 


be expected that k will be larger for a continually agitated, well-mixed bed as opposed to a quiescent, 
stationary bed. For a well-mixed bed, it may be expected that the evaporation of the condensed water 
occurs on the surfaces of the individual particles of the regolith so that S in this case is proportional to the 
surface area of the regolith particles. For a given amount of regolith, the surface area will be larger for 
smaller particles. Thus, there may be some benefit to load the reactor with relatively smaller particles. For 
the case where the regolith bed is stationary the effective surface area is expected to be related to the 
cross-sectional area of the reactor. In this case, evaporation may take place starting from the surface of the 
bed and progressing through the bed in layers.The spatial variation in condensed water fraction then has 
to be taken into account which will modify the analysis presented in the previous section. 

For the case where the reactor temperature increases with time, Equation (10) must in general be 
solved numerically. Even if the parameters V, k, and S are taken as constant during the heating, the 
saturation vapor density, p„ remains a strong function of the temperature. In order to assess the effects of 
the time constant for the case of varying temperature, Equation (10) has been solved for the case of a 
linear increase in temperature with time. The corresponding variation of reactor pressure is plotted in 
Figure 2 for two different values of the time constant (i.e., t, = 100 s and t, = 500 s, respectively). The 
temperature is assumed to increase linearly from 298 K (25 °C) to 423 K (150 °C) over 1 hr. The vapor 
pressure-temperature relationship needed for the calculations as expressed by Equation (1) have been 
modeled using the Hyland-Wexler correlation (Ref. 10) for liquid water. 

As in the case of constant temperature, it is found that decreasing the time constant results in a more 
rapid rise in pressure. Comparing Figure 1 and Figure 2, it is also seen that for the case of varying 
temperature the pressure -time curve is concave upward whereas for the case of constant temperature the 
curve is concave downward. This shows that the rate of increase of vapor pressure in the reactor increases 
as the temperature rises whereas for the constant temperature case, the rate decreases with time as the 
vapor pressure builds inside the reactor. 

The behavior of the vapor pressure -time curve predicted by the model and depicted in Figure 2 has 
been observed in experiments conducted at the NASA Johnson Space Center (JSC). In these experiments, 
JSC-1A, a lunar regolith simulant, with added water was placed in a small covered crucible and heated. 
The increase in pressure in the crucible was monitored. An example of the observed pressure-time 
behavior is shown in Figure 3. The experimental curve is seen to exhibit the concave up behavior similar 
to the model prediction. During the observation time the crucible temperature increased from 
approximately 273 to 423 K. Further quantitative comparisons with the limited crucible experiments have 
not been carried out as the experimental efforts have now shifted to the development of the tubular reactor 
configuration for the RESOLVE experiment. 
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Figure 3. — Experimentally observed increase in vapor pressure with 
time during heating. 



Figure 4. — Nondimensional water evaporation rate as a function of 
temperature for two different effective density parameter values. 

Effects of Reactor Temperature and Pressure on Evaporation During Venting 

According to Equation (1), the saturation pressure P s is a nonlinear function of the reactor temperature 
T r . Assuming that the free volume of the reactor and the evaporation time constant do not change during the 
vent, Equation (16) indicates that the water evaporation rate is also a nonlinear function of the reactor 
temperature. This is shown in Figure 4 which plots the nondimensional evaporative mass rate as a function 
of the reactor temperature for two different values of the effective density parameter (m 0 -mi)/V. The 
corresponding values of the nondimensional evaporation time are plotted in Figure 5. Increase in the 
effective density parameter indicates a smaller value of the free volume for a given amount of condensed 
water in the reactor. Equation (18) shows that the evaporative mass flux rate is then reduced. 

The water evaporation rate is proportional to the difference between the saturation vapor pressure, P s , at 
the reactor temperature T r , and the regulated reactor pressure P r . If P, is much less than P s then the regulated 
pressure value is not important in the evaporation process. This is shown in Figure 6 which plots the 
nondimensional evaporation time as a function of reactor temperature for two different values of the 
regulated reactor pressure, i.e., 0.01 and 0. 1 atm, respectively. The density parameter for these calculations 
is taken to be 100 kg/m 3 . At the lower temperatures, the evaporation time is relatively smaller for the 0.01 
atm case. However, beyond approximately 350 K as the saturation pressure becomes much larger than 0.1 
atm, there is an insignificant difference between the evaporation times for the 0.01 and 0. 1 atm cases. 
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Figure 5. — Nondimensional water evaporation time as a function of 
temperature for two different effective density parameter values. 



Figure 6. — Nondimensional water evaporation time as a function of 
temperature for two different regulated reactor pressure values. 


Conclusion 

A model has been developed to analyze the processes occurring during the heat-up of icy regolith and 
extraction of the evolved water vapor. The first part of the model calculates the vapor pressure increase in 
a closed cell containing the initially wet regolith as the temperature is raised and the second part analyzes 
the evaporation time and mass evaporation rate during reactor venting. The flux of water vapor from the 
evaporating surface to vacuum is first obtained and then the external vapor pressure is accounted for. The 
rate of increase of the vapor pressure during heat-up depends upon a characteristic time constant related to 
the free volume in the cell, the effective surface area for vaporization and a mass transfer coefficient. 
Reducing the time constant increases the rate of vaporization. The shape of the vapor pressure-time curve 
obtained from the model is shown to be similar to that obtained experimentally. The evaporation time 
during the vent period depends strongly on the regolith temperature but is not sensitive to the pressure to 
which the reactor is vented as long as this pressure is much less than the saturation vapor pressure at the 
regolith temperature. 
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